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Abstract 
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threshold, in a Bethe-Salpeter model. We make predictions for the radiative decay 
widths of the B c excited states. We compare our results with those of other models. 
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1 INTRODUCTION 



The B c meson discovered by the CDF collaboration [1] in pp collisions at y/s = 1.8 TeV 
completes the family of mixed flavor mesons. The B c meson has a b anti-quark and a c 
quark. Current and future experiments at the Tevatron and LHC are expected to provide 
large samples of the excited states of the B c mesons [2]. This will make possible the study 
of the spectroscopy and the decays of the B c mesons. The B c meson family lies intermediate 
in mass and size between the cc (J/ip) and the bb (T) families where the heavy quark 
interactions are believed to be understood rather well. Comparison between experimental 
measurement and theoretical results will improve our understanding of these interactions 
and guide us in the search for multiquark and molecular exotics such as the recently claimed 
(discovered) D sJ [3, 4, 5] and X(3872) [6]. 

Different models [7, 8, 9, 10, 11, 12] including various versions of potential models and 
QCD sum rules have been used to evaluate the B c spectrum yielding results consistent 
with the experimentally measured ground state mass and lifetime. The B c mesons have 
non-vanishing flavor quantum numbers which are conserved in strong and electromagnetic 
interactions. Therefore, the B c states, below the open flavor BD threshold, can only decay 
weakly or radiatively. These states are expected to be relatively long-lived and easier to 
be observed experimentally. From the theoretical side, weak and radiative decays are free 
from uncertainties encountered in strong decays which makes the decays of these states 
theoretically more tractable. 

In a previous paper [13] , we have evaluated a limited set of the B c spectrum using a model 
based on reductions of the Bethe-Salpeter equation (BSE). We have used a set of parameters 
fixed from previous investigations of other meson spectra. Our results agreed very well with 
the experimentally measured ground state mass and lifetime. We also evaluated the B c decay 
constant, the b antiquark and the c quark inclusive decay widths and the weak annihilation 
width. 

We also evaluated the exclusive semileptonic (B c — > P{V)ev) and two-body nonleptonic 
(B c — > PP, PV, VV) decay widths [14], where P (V) denotes a pseudoscalar (vector) meson. 
We used the BSE amplitudes to evaluate the semileptonic form factors and used factorization 
to obtain the nonleptonic decay widths in terms of the semileptonic form factors and the 
weak decay constants. 

In the present paper, we evaluate the complete B c spectrum below the open flavor BD 
threshold and consider the radiative El and Ml electromagnetic transitions. This comple- 
ments our picture [13, 14] of the B c mesons. Radiative decays are the dominant decay modes 
of the B c excited states having widths of about a fraction of MeV, much greater than the 
weak widths at the order of meV. Therefore, accurate determination of the masses and the 
radiative decay widths will be extremely important for understanding the B c spectrum and 
distinguishing exotic states. 

The paper is organized as follows. In the next section we briefly outline our model 
and compare our spectrum with those of other models. We then evaluate the El and Ml 
radiative decays. Finally we discuss our results. 
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2 MODEL AND SPECTROSCOPY 



We applied a relativistic model based on reductions of the BSE to evaluate the B c spectrum. 
The BSE is a suitable starting point for treating hadrons as relativistic bound states of 
quarks and antiquarks, just as the Dirac equation provides a relativistic description of a 
fermion in an external field. The BSE for a bound state may be written in momentum space 
in the form [15] 



g-\p,pW(p, p ) = J ^±^v(P, P -p')4>(p, P ')d 4 P ' (l) 

Where P is the four-momentum of the bound state, p is the relative four-momentum 
of the constituents. The BSE has three elements, the two particle propagator (G) and the 
interaction kernel (V) which we provide as input, and the amplitude ijp) obtained by solving 
the equation. We also solve for the energy, which is contained in the propagator. We used 
a reduction of the BSE where the two particle propagator is modified in a way that keeps 
covariance and reduces the four-dimensional BSE into a three-dimensional equation [16]. 
We considered an interactional kernel that consists of two terms, one for the short range 
one gluon exchange Voge and the other for the long range phenomeno logical confinement 
interaction Vcon [17]. 



Voge + V CO N = ~^ s - — + ahm — — — — -. 2 

3 ip — p) ^oo/i z —{p — p') z + n 2 

Here, a s is the strong coupling, which is weighted by the meson color factor of |, and 
the string tension a is the strength of the confining part of the interaction. While the one 
gluon exchange Voge has the vector nature, we adopt a scalar Lorentz structure for Vcon 
as discussed in [16]. We solve for the energies and the amplitudes in momentum space and 
transform these amplitudes into coordinate space. 

We have included seven parameters in our model, four masses (m u = morris, m c ,m^}, 
two parameters to fix the strong coupling ct s and control its running with the meson mass, 
and the last parameter is the string tension a of the confining interaction. We fixed the 
parameters of our model by fitting the spectra of other mesons as described in [17]. We 
obtained a good fit for a wide range of meson masses with root mean square deviation from 
experimental masses of about 50 MeV. 

Table 1 compares the parameters relevant to the B c mesons of our model with those of 
different models in the literature. In Table 1, m c and rrib are the masses of the c and b 
quark respectively, while a s is the strong coupling of the one gluon exchange and a is the 
string tension of the confining interaction. In many models, including ours, a s runs with the 
meson mass, thus Table 1 gives a s at the scale of the ground state mass of the B c mesons. 
The model used in [8] employs the Martin potential [18] which is not linear but varies with 
powers of the quark antiquark distance. 
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We notice that our m c and m& values are smaller that those of other models, while our 
string tension a is larger. The values of the strong coupling a s are consistent around 0.36 
except in [11] where a s is 0.265 and in [12] where a s is 0.21. 



Table 1: The parameters of different models relevant to the B c mesons. 





This work 


EQ [7] 


GKLT [8] 


EFG [11] 


GI [12] 


m c (GeV) 


1.39 


1.48 


1.8 


1.55 


1.628 


m b (GeV) 


4.68 


4.88 


5.174 


4.88 


4.977 


a s 


0.357 


0.361 


0.391 


0.265 


0.21 


a (GeV 2 ) 


0.211 


0.16 




0.18 


0.18 



Fig. 1 shows the B c spectrum of our model. The horizontal dashed line represents the 
BD threshold (7143 MeV). States above the BD threshold can decay strongly into two 
heavy-light mesons while those below that line can only decay weakly or radiatively. Our 
result for the B c ground state mass (6.380 GeV/c 2 ) agrees very well with the Experimental 
result of the CDF collaboration 6.40 ± 0.39 (stat.) ± 0.13 (syst.) GeV/c 2 [1]. 
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;ure 1: 


The B c 
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Table 2 compares our B c spectrum with those of some other models. One may notice 
that the hyperfine splitting of our model ( 3 S\ - 1 Sq difference) is smaller than those of 
other models, while the fine splitting of the P states is larger in our model. We treat the 
interactions responsible for these splittings directly in the bound state problem while Eichten 
and Quigg [7], for example, treat them perturbatively. Experimental results for the excited 
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states of the B c mesons are needed to clarify these differences and improve our knowledge 
of the parameters shown in Table I. 

Table 2: B c spectrum in units of GeV. 
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1 3 £>2 


7.000 


7.012 


7.001 


7.077 


7.041 


1^2 


7.001 


7.009 


7.016 


7.079 


7.036 


1 3 ^3 


7.003 


7.005 


7.007 


7.081 


7.045 


2 3 P 


7.081 


7.108 


7.088 


7.091 


7.122 


2 3 P X 


7.136 


7.135 


7.113 


7.126 


7.145 


2^! 


7.139 


7.142 


7.124 


7.145 


7.150 


2 3 P 2 


7.186 


7.153 


7.134 


7.156 


7.164 


3% 


7.198 


7.244 




7.193 


7.250 




7.215 


7.280 




7.235 


7.272 


4% 


7.452 


7.562 






7.572 


4 3 S' 1 


7.468 


7.594 






7.588 



We note that the lifetime of the B c ground state as reported by CDF is 0.46^o;Jg (stat.) 
± 0.03 (syst.) ps [1] corresponding to a width of about 1.43 meV. Three different processes 
contribute to the B c ground state width : inclusive b decays, inclusive c decays and b-c 
annihilation. For excited B c states the corresponding weak decays will have partial widths 
similar to or less than the ground state widths (since the weak b-c annihilation will be 
forbidden for S^O states). Radiative decays having widths of the order of a fraction of MeV 
(as evaluated in the next section) will be the dominant decay modes of the B c excited states. 

Using the B c spectrum and the BSE amplitudes in momentum space or the transformed 
coordinate space amplitudes one can evaluate the radiative transition widths between the 
B c states. This is what we address in the next section. 
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3 RADIATIVE DECAYS 



The electromagnetic radiative decays of the B c mesons are of the electric dipole (El) and 
the magnetic dipole (Ml) types. The El partial decay widths can be written as [7, 19, 20] 

r*i(i - / + 7) = 4Q< 3 eq>2 (2J / + l)u 3 \{f\r\i)\ 2 C fi (3) 
where the mean charge is 

< m b e c -m c e- b 
m h + m c 

e c = 2/3 is the c quark charge and ei = 1/3 is the charge of the b antiquark in units of |e|, 
a is the electromagnetic fine structure constant, u is the photon energy 

U = ^m (5) 

and the statistical factor C/j is given by (S = Si = Sf) 

C fi = max(L i , L f ) {^^f . (6) 

The mean charge expresses the fact that the emitted 7 can be attached to the c quark or the 
b antiquark. The lighter c quark is more efficient in this process. The statistical factor C/j 
results from the angular momentum coupling. Table 3 compares our results (This work) and 
the results of Eichten and Quigg (EQ) [7] for the transition energy (in MeV), the transition 
matrix element |(/|r|i)| (in GeV" 1 ), and the transition width (in keV). We notice that the 
transition energies may differ by a factor of about two while the transition matrix element 
are very close to each other. The transition matrix elements in our model depends on the 
initial and final values of L, J, S since our model treats the spin-orbit, the spin-spin, and 
the tensor interactions in the bound state problem while the transition matrix elements of 
Eichten and Quigg [7] depend on the initial and final L only since they take care of these 
splittings perturbatively. 

The Ml partial decay widths between the S-wave states can be written as [7, 19, 20] 

rMi(^/ + 7) = ^(2J/ + lV 3 |(/|jo(W2)K)| 2 , (7) 
where the magnetic dipole moment is 

m^ec — m c e-^ 
Am c mb 

Allowed Ml transitions correspond to triplet-singlet transitions between S'-wave states of 
the same n quantum number, while hindered Ml transitions are either triplet-singlet or 
singlet-triplet transitions between S'-wave states of different n quantum numbers. In the non- 
relativistic limit these hindered transitions are suppressed by wave function orthogonality. 
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Table 3: Comparison of the results for the El transition rates of the B c mesons. The columns 
labelled (This work) present our results. 





oo (MeV) (GeV" 1 ) 


r(<- 


->/ + 


7) (keV) 


Transition 


This work EQ[7] This work EQ[7] 


This 


work 


EQ[7] 



l 3 P 2 ^l 3 ^i+7 408 397 1.613 1.714 109.8 112.6 

l 3 P 1 ^l 3 ,S 1 +7 347 382 1.619 1.714 67.8 99.5 

l 3 P 1 ^l 1 S + 7 381 450 1.537 1.714 0.0 0.0 

l 1 P 1 ^l 3 S 1 +7 349 387 1.615 1.714 0.0 0.1 

-> l 1 ^ + 7 383 455 1.531 1.714 81.8 56.4 

l 3 P ^l 3 ^i+7 270 353 1.617 1.714 32.1 79.2 

2 3 Si^l 3 P 2 + 7 58 151 1.870 2.247 0.7 17.7 

2 3 Si^l 3 Pi+7 123 167 1.836 2.247 3.9 14.5 

2 3 Si^l 1 P 1 +7 120 161 1.887 2.247 0.0 0.0 

2 3 Si -> 1 3 P + 7 200 196 1.862 2.247 5.8 7.8 

2 1 S'o ^ l 3 Pi +7 102 125 1.915 2.247 0.0 0.0 

2 1 S , ^l 1 Pi+7 99 119 1.965 2.247 7.1 5.2 

i 3 P> 3 ^i 3 P 2 + 7 163 258 2.404 2.805 18.7 98.7 

i 3 P> 2 ^l 3 P 2 + 7 160 258 2.405 2.805 4.4 24.7 

1 3 £> 2 -> I 3 Pi + 7 224 274 2.395 2.805 35.8 88.8 

l 3 P> 2 ^l 1 Pi+7 221 268 2.436 2.805 0.0 0.1 

l 3 P> 1 ^l 3 P 2 + 7 121 258 2.431 2.805 0.2 2.7 

l 3 P> 1 ^l 3 P 1 +7 184 274 2.415 2.805 11.4 49.3 

l 3 P> 1 ^l 1 P 1 +7 182 268 2.454 2.805 0.0 0.0 

1 3 P>! -> 1 3 P + 7 262 302 2.434 2.805 43.9 88.6 

l 1 P> 2 ^l 1 Pi+7 222 268 2.433 2.805 120.8 92.5 

2 3 P 2 ^l 3 1 Si+7 729 770 0.194 0.304 9.1 25.8 

2 3 P 2 ^2 3 S 1 +7 285 249 2.525 2.792 91.3 73.8 

2 3 P 2 ^l 3 P> 3 + 7 181 142 2.249 2.455 31.5 17.8 

2 3 P 2 -> 1 3 P 2 + 7 184 142 2.238 2.455 5.8 3.2 

2 3 P 2 -> l 3 Pi + 7 223 142 2.046 2.455 0.6 0.2 

2 3 P 1 ^l 3 ,Si+7 684 754 0.178 0.304 6.3 22.1 

2 3 Pi -> 2 3 Si + 7 236 232 2.553 2.792 53.5 54.3 

2 3 P!^l 3 P 2 + 7 136 125 2.277 2.455 12.0 9.8 

2 3 P! -> l 3 Dx + 7 175 125 2.078 2.455 7.2 0.3 

2 1 Pi^l 3 Si+7 686 760 0.209 0.304 0.0 2.1 

2 1 P 1 ^2 3 Si+7 239 239 2.509 2.792 0.0 5.4 

2 1 Pi^l 3 P> 2 + 7 138 131 2.232 2.455 0.0 11.5 

2 1 P 1 ^l 3 P> 1 +7 177 131 2.029 2.455 0.0 0.4 

2 3 P ^l 3 1 S 1 +7 634 729 0.193 0.304 5.9 21.9 

2 3 P ^2 3 1 Si+7 183 205 2.531 2.792 24.3 41.2 

2 3 P ^l 3 P>i+7 121 98 2.053 2.455 9.3 6.9 
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Taking relativistic effects in the wavefunction (in addition to the large a; 3 dependence of 
the width) makes the widths of these transitions comparable with the widths of the allowed 
ones. Table 4 compares our results (This work) and the results of Eichten and Quigg (EQ) 
[7] for the transition energy (in MeV), the transition matrix element |(/|jo(^"/2)|i)|, and 
the transition width (in keV). 

Table 4: Comparison of the results for the El transition rates of the S'-wave B c mesons. The 
columns labelled (This work) present our results. 







lo (MeV) 


\(f\jo(kr/2)\i)\ 


r(i_/ + 7 ) (keV) 


Transition 




This work 


EQ[7] 


This work 


EQ[7] 


This work EQ[7] 


2 3 S 1 -> 2% 


+ 7 


21 


43 


0.9949 


0.9990 


0.0037 0.0289 


2 3 S 1 -> l 1 ^ 


+ 7 


496 


606 


0.0523 


0.0395 


0.1357 0.1234 


2^5*0 — ^ 1 3 jSi 


+ 7 


443 


499 


0.0393 


0.0265 


0.1638 0.0933 


l 3 ^ - l 1 ^ 


+ 7 


36 


72 


0.9979 


0.9993 


0.0189 0.1345 



Table 4 shows differences between the transition energies (spin-spin splittings), while the 
transition matrix elements are comparable. It also shows that hindered transitions have 
widths at the same level as the allowed ones. 

4 SUMMARY 

We have evaluated the B c spectrum below the BD threshold using a reduction of the BSE. 
We have made predictions for the transition rates of the El and Ml radiative decays. We 
compared our results with the results of other models in the literature. Experimental results 
will help clarify the spin-spin and spin-orbit splittings of different models and consequently 
improve our knowledge of physical quantities such as the quark masses and the strong cou- 
pling a s at the scale of Mb c - Measurements of radiative transitions and comparison with 
such results may indicate the existence of exotic multiquark or molecular exotics. 
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